Cyclic phase in F=2 spinor condensate: Long-range order, kinks, and roughening 

transition 
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We study the effect of thermal fluctuations on homogeneous inflnite Bose-Einstein condensate 
with spin _F = 2 in the cyclic state, when atoms occupy three hyperfine states with mp = 0, ±2. 
We use both the approach of small-amplitude oscillations and mapping of our model on the sine- 
Gordon model. We show that thermal fluctuations lead to the existence of the rough phase in one- 
and two-dimensional systems, when presence of kinks is favorable. The structure and energy of a 
single kink are found. We also discuss the effect of thermal fluctuations on spin degrees of freedom 
in F = 1 condensate. 

PACS numbers: 03.75.Kk, 03.75.Lm 



I. INTRODUCTION 

Recent progress in the experimental physics of cold 
atoms allows for the creating of quasi-two-dimensional 
atomic gas either using one-dimensional optical lattices 
or applying tight axial trapping [ll Q, S 13 ■ Many prop- 
erties of low-dimensional systems differ from that for 
the three dimensional ones. For instance, it is well- 
known that Bose-Einstein condensation is impossible in 
two-dimensional {2D) situation at any nonzero tempera- 
ture, since long wavelength thermal fluctuations destroy 
a phase coherence. However, at low temperatures, there 
is still a quasi-long-range order, which disappears at the 
temperature of Berezinskii-Kosterlitz-Thouless (BKT) 
transition. Low temperature phase is often called qua- 
sicondensate, and this conception was introduced to the 
context of physics of alkah atoms gases in Refs. [l,|a|, see 
also Ref. ,71 ■ In was shown in those works that, in finite 
low-dimensional systems, there is still a long range order 
at low temperatures, and one has a true condensate, but 
at higher temperatures, coherence is again lost, and we 
have a quasicondensate. 

The aim of the present paper is to study thermal fluc- 
tuations in homogeneous spinor condensates. We are 
mainly interested in the cyclic phase of i^ = 2 conden- 
sate, since it has unusual properties like phase locking 
and kinks, which are absent in other phases of i^ = 2 
condensate. In the recent experiments H, y, llfl HH Hjl j 
F = 2 spinor Bose-Einstein condensates have been cre- 
ated and studied. The order parameter in F = 2 system 
has five components; and, in the cyclic phase, all the par- 
ticles populate three hyperfine states with mp — 0, ±2. 
The characteristic feature of the cyclic phase is the fact 
that the ground state energy depends on the relative an- 
gle X among the phases Sj of different components of the 
order parameter, x = 2So — S2 — S-2, through the spin- 
mixing term. This leads to the peculiar phase locking 
phenomena, since the energy has a term proportional to 
cos X- In this paper we concentrate on thermal fluctua- 
tions of phases Sj of different components of the order 
parameter as well as of x- First, we solve the Bogoliubov- 
de Gennes equations and calculate the mean square fluc- 



tuations of Sj. Each individual phase Sj in this case 
behaves in the way similar to the case of scalar system. 
Namely, there is no long-range order for Sj in one and 
two dimensions, whereas this order is kept for 3D sit- 
uation. At the same time, the long-range order is still 
preserved for x- After this, we note that the Bogoliubov- 
de Gennes equations can be insufficient for the analysis 
of thermal fluctuations in the cyclic state, since this ap- 
proach uses an expansion of the energy in the vicinity 
of one of the infinite number of equivalent minima at 
X = 2-kI, I being an integer number, and does not take 
into account a global periodic structure of the energy 
in functional space. In other words, cosx contribution 
to the energy is changed by the quadratic potential well 
with the infinite height. One of the possible solutions of 
the Gross-Pitaevskii equations for the cyclic phase is a 
kink, which separates two spatial domains with x differ- 
ent by 27r from each other. At zero temperature, kink is 
energetically unfavorable. However, at finite temperature 
formation of kinks can become favorable due to entropic 
reasons. We determine analytically the structure and en- 
ergy of a single kink. Also we note that our model can 
be mapped on the well-known sine-Gordon model. For 
this model, it was established before that, in low dimen- 
sions, a so called roughening transition can occur, when 
the system becomes unlocked from one of the minima of 
X- In ID situation, a long range order for x is absent for 
any nonzero temperature due to thermally-excited kinks, 
and the system is rough. If the system has a finite length, 
a nonzero temperature appears, below which still there is 
a long range order for x- We estimate this temperature 
from simple entropic arguments. In two dimensions, a 
roughening transition occurs at some finite nonzero tem- 
perature. From known results for the sine-Gordon model 
and some simple qualitative considerations we show that 
the temperature of the roughening transition is close to 
the BKT critical temperature and possibly coincides with 
it. We also discuss briefiy the case of ferromagnetic F = 1 
condensate. Based on the solutions of the Bogoliubov-de 
Gennes equations, we find that there is no long-range 
order in the one- and two-dimensional systems for the 
direction of spin. In finite systems, the order still exists 



at low temperatures, similarly to the phase coherence in 
scalar condensates. 

The paper is organized as follows. In Section II we 
solve the Bogoliubov-de Gennes equations for _F = 2 
cyclic phase and find correlation functions. In Section 
III we determine the structure and energy of a single 
kink. In Section IV we discuss the order-disorder tran- 
sition due to the proliferation of kinks. In Section V we 
study the case of ferromagnetic F = I condensate. We 
conclude in Section VI. 



II. CYCLIC PHASE: SMALL AMPLITUDE 
OSCILLATIONS 

We consider infinite homogeneous F = 2 condensate 
with a given density of particles n in zero magnetic field. 
The energy of the system depends on three interaction 
parameters a, /?, and 7, which can be defined as |13Lll4| 



y(4g2 + 3.g4), 



/3 = -y(52 



54), 



7= g(.9o-54)- ^(32-54), 



where {q = 0,2,4) 



AttH^ 



9q 



(1) 
(2) 
(3) 

(4) 



and aq is the scattering lengths characterizing collisions 
between atoms with the total spin 0, 2, and 4. In real 
atomic condensates, a^ 13, 7. 

The order parameter in i^ = 2 case has five components 
'^i ii = —2, —1, 0, 1, 2). The energy of the system is given 
by [Hill 



F ^ dv 



a 



Tr-^*^^j + ^^^n^^j^k 



2m 






(5) 



where integration is performed over the system volume, 
repeated indices are summed. Fa (a = x,y,z) is the an- 
gular momentum operator, which can be expressed in a 
usual matrix form. 

In the absence of magnetic field and rotation, conden- 
sate can be in three different states [13, as seen from 
Eq. (5). These states are called ferromagnetic, cyclic and 

polar I13. In cychc phase, *±i — and *_2 



Vn A0 



where 9 is an arbitrary phase 



(energy of the system is degenerate with respect to 9). 
The ground state value of x in these notations is zero. 



Depending on values of scattering lengths aq ferror nag - 
netic, cychc or polar phase has the lowest energy [l3| . 
Extended Gross-Pitaevskii equations can be obtained as 
usually from the condition of minimum of free energy of 
the system Eq. (5). In the cyclic state, the last term 
in the right-hand side of Eq. (5) contains a spin-mixing 
contribution: 



i^sm=7l*0r|*2||^-2|cOSX. 



(6) 



Now we analyze small-amplitude oscillations of the or- 
der parameter in the cyclic phase. The fluctuations of 
phase of the given component of the order parameter 
can be expressed through the fluctuations of the order 
parameter itself as 6Sj = Im{5'^ j) / \^ j\. For 5x, we 
have 5x — 2SSq — SS2 — SS^2- 

The deviations of the five components of the order pa- 
rameter from their equilibrium values (5^j can be repre- 
sented as 

S'^j = y^ \c^^^ exp(ikr -I- ioji^t) 

k 

+d\^''* exp{-ikr - ioji^t)^ , (7) 

(i) (i) 

where c^ and d^ are constants to be found. Fluctu- 
ations of ^±1 are decoupled from that for ^0, ^±2, as 
follows for the Bogoliubov-de Gennes equations for the 
cyclic phase. The spectrum for oscillations of ^0, ^±2 
has three branches. For the first time, it was obtained in 
Ref. [13, and our results coincide with those results: 



huj, 



(1) 



2u2 



2m 



2-fn, 



(2) 



hu), 



^f = 



V 2m J 



2an- 



2m '' 



\ 2m ) 



■4:(3n 



fi2k2 

2m ' 



(8) 



(9) 



(10) 



The thermal distribution of quasiparticles is given by 
(j - 1, 2, 3) 



N, 



U) 



1 



eMf^^^^lkBT)-! 



(11) 



One can also find eigenvectors corresponding to the 
branches (8)-(10). Eigenvectors for Eqs. (9) and (10) 
give no contribution to (5%, as it can be seen from the 
solutions of the Bogoliubov-de Gennes equations. This is 
the first branch (8), which is responsible for the fluctu- 
ations of X- Physically, this is due to the fact that only 
branch (8) depends on the value of 7, which yields the 
potential well for x via the spin-mixing term (6). Finally, 
we obtain simple expression for 5x'- 



6x 



1 



y {exp(ikr + iuj\st) 



2i^/V^ ^ 

— exp(— ikr — zwki)} 



(12) 



Next we can determine the behavior of the mean square 
fluctuations of x at large distances. We found that this 
quantity in all dimensions tends to a constant, when the 
distance tends to infinity. This means that there is a 
long-range order for x. The reason is the presence of 
the gap in the excitation energy (8). As an example, 
we present below the derivation and results for the 2D 
situation. From Eq. (12) we have 



{5x{r) - -^xlO))' = - -^ 5] Ak, Ak, , (13) 



Wn 



ki,k2 



where 



Av 



exp(icj|j. 't) (exp(ikr) — 1) 

— exp(— JWi^ h) (exp(— ikr) — 1) , 



After averaging over the time, only terms with fci = ^2 
survive in the expansion (13) and for the correlator we 
have 



J 



[5x{v) - 5xm 



E ^. 



0) 






^ikr cos ipi 



^ikr cos :p2 



im 



where ipi and ip2 are angles between r and ki and k2, 
respectively. Now we can switch from the integration to 
the summation, and after the integration over (pi and ip2 
we get 



([6x{r) ~ Sx{0)f)^ - ^ fkdkNi''' (1 - Jo(fcr)) 



(15) 
where Jo{r) is the Bessel function. This integral can be 
evaluated analytically at large r and in the limits of low 
and high temperatures. At low temperatures, kgT <C 



in, 



N, 



(1) _ 



exp(- 

(1) 



keT > 7n, N^'> = keT/nw)^ 



hLul^^/ksT). 
,0) 



At high temperatures. 
Finally, we have 



[dx{r) ~ 6xm 



T/TaeM-^iri/kBT), 
T/Td\n{kBT/2-fn), 



ksT < 7n, 
ksT > 7n, 



(16) 



r 



where Td = 2iTh'^n/kBm is the temperature of quan- 
tum degeneracy. 

The second and third branches of the spectrum, given 
by Eqs. (9) and (10), are responsible for the fluctuations 
of phases of individual components of the order param- 
eter, 5*0, S±2. Our calculations of the mean square fluc- 
tuations of these phases revealed the behavior, which is 
similar to the scalar condensate. The long-range order is 
absent in one- and two-dimensional situations. It is in- 
teresting to note that, at the same time, the long-range 
order exists for the linear combination of these phases, x- 
Fluctuations of populations of different hypcrfine states 
are small at T ^ T^. 



III. KINK: STRUCTURE AND ENERGY 



lution to the extended Gross-Pitaevskii equations. How- 
ever, these equations have also a nonhomogeneous solu- 
tion corresponding to the kink or domain wall between 
two spatial regions with values of x different by 271"/ from 
each other, since the energy of the system is degenerate 
with respect to I. Kinks are known to play an important 
role in the physics of low-dimensional systems, see e.g. 
the textbook Te| . Now we find the structure and energy 
of the one-dimensional kink in the cyclic phase for the 
most important case, when I — 1. We assume that all 
the quantities depend only on one coordinate x. For x 
we have a boundary conditions x(c«) = ^tt, x(— oo) — 0. 
The phase-dependent part of the energy consists of two 
contributions, one is the kinetic energy, and another one 
is the spin- mixing term (6): 



In the previous Section we have studied small- 
amplitude oscillations in the vicinity of homogeneous so- 



J 



F, 



ph 



-t-oo 



dx 



2m 



J2 l*2jl' (^^2,)' + 7 l*o|' |*2| |*-2| cosx 



i=-i 



(17) 



We can rewrite this expression in the diagonal form as 



/' + 00 



F, 



ph 



dx 



h^nl 



jn 



n lb -L 2 7" 



h'^n / 1 ,„ ^ 



r 



-V5. 



-2? + -^ (V52 



- vs^2 + zys-o) 



(18) 



One can easily see from Eq. (18) that the minimum 
energy solution for the kink corresponds to the conditions 
VS2 - VS-2 = and VS'a + VS-2 + '^.V So = 0. In 
this case, the energy (18) can be mapped on the sine- 
Gordon energy functional. The function x{^) is found 
from the condition of minimum of Fph'- 5Fph/Sx = and 
boundary conditions. Finally, we have 



S2{X) 



= 27r-4tan"^exp(-a;/^x), (19) 

= S.2{x)=~So{x)^~jx{x), (20) 



where 



I h^ 1 

2m 4771 



(21) 



is the healing length for x- It also gives a characteris- 
tic length of the kink. The spatial structure of a single 
kink is presented in Fig. 1, where we have plotted x- 
dependences of Xj 'S'±2, and Sq. A single kink moving in 
the space is a soliton. As seen from Eq. (18), the energy 
of the kink Fkink scales as 7^;,. ~ ^. More accurate 
calculation based on Eq. (19) yields 



kink 



,3/2 



^27 



(22) 



Note that in real atomic condensates 7 has to be much 
smaller than a and, therefore, ^^^ should far exceed the 
coherence length, which has a meaning of a length-scale 
for the density modulations. 

The same structure of the kink can be obtained from 
the extended Gross-Pitaevskii equations supplemented 
by the boundary conditions for x- These equations are 
rather cumbersome and we do not present them here (see 
e.g. our previous work [l^)- It is easy to see by the direct 
substitution that Eqs. (19)-(21) with the constant popu- 
lations of each magnetic sublevel yield the exact solution 
to these equations. 



IV. ROUGHENING 

At zero temperature, the sine-Gordon system in any 
dimension is locked in one of equivalent minima corre- 
sponding to X = 2ttI. At nonzero temperature, presence 
of kinks can be favorable due to their entropic contribu- 
tion to the free energy [I3] ■ Kinks destroy the long-range 
order and lead to the rough phase, in which height of 



fluctuations diverges with tending the system size to in- 
finity. In one dimensional system of length L (L ^ ^^), 
the entropy corresponding to the single kink can be esti- 
mated as InL/^^, where L/(^^ is just a number of places 
to put a kink. Therefore, a roughening temperature is 



rplD 



F, 



kink 



ks In L/^x 



(23) 



In the case of infinite ID system (L -^ 00), T^ = 0; 
and this result is very different from the result of small- 
amplitude oscillations approach, which shows a presence 
of long-range order for x- Relation (23) is well-known for 
the sine-Gordon model and can be derived using more rig- 
orous analysis, see e.g. Ref. [^^ and references therein. 
Note that in ID finite system, T^^ ^ 0, as 7 ^ 0. 

In two dimensions, a roughening temperature T^^ can 
be also estimated in a simple manner. In this case, Fkink 
given by Eq. (22) is the energy of a unit length of a kink. 
The 'critical nucleus' for the formation of the rough phase 
is represented by the step for x of height 27r in a func- 
tional space with a perimeter « 27r^^ . The total energy of 
such a nucleus is of the order of 2TrFkink£,x — T^fi^ns/2m, 
where we have changed density of atoms n by super- 
fluid density Us just below the roughening transition. 
It is interesting to note that this energy is indepen- 
dent on 7. A roughening temperature T^ can be es- 
timated by equating the nucleus energy to the thermal 



one. 



Tr 



Trh^ns/2kB'm. A more rigorous calculation 
for the 2D sine-Gordon model, based on renormalization 
group analysis [13], yields the same result: 



T, 



2D 



2fcBTO' 



(24) 



which is independent on 7. The expression (24) for the 
roughening temperature coincides with the well-known 
result for the temperature of the BKT transition, see 
e.g. |3- To describe properly the behavior of the system 
in this strongly fluctuative region at high temperatures, 
one needs more careful and detailed analysis. However, 
from the considerations presented here, we can conclude 
that a long range order for x in 2D systems survives up to 
quite high temperatures (of the order of the BKT critical 
temperature) . Possibly, the roughening transition occurs 
simultaneously with the BKT transition, and both types 
of topological defects, kinks and vortices, proliferate to- 
gether. 

If the size of ID or 2D finite system is less than f^ 
(which can be much larger than the coherence length), 
a formation of a kink is impossible, and, therefore, we 



expect that x is nearly constant inside the cloud. At 
the same time, a value of x can be different from 27:1 
due to thermal fluctuations. We have studied the similar 
situation before in Ref. zl] , where the case of harmon- 
ically trapped quasi-two-dimensional F — 1 condensate 
was treated. An infinite homogeneous F = 1 condensate 
in zero magnetic field can be either in polar or ferro- 
magnetic states, where atoms populate only one or two 
hyperfine states and the spin-mixing terms in the energy 
are zero in the equilibrium. However, in the case of a 
trapped rotated condensate containing vortices, there are 
some regions on the phase diagram, where all the three 
hyperfine states are populated, and the ener gy dep ends 
on the relative phase among them pi I2II2I 1211261 l27t . 
This vortex phases can be both locally and globally sta- 
ble, and some of them, like Mermin-Ho vortex, have an 
axial symmetry, which makes them rather simple. An- 
other method to create such states was recently used ex- 
perimentally in Ref. |2q , where a microwave energy was 
injected to the system leading to the redistribution of the 
particles from spin —1 state to spin and 1 states. We 
see that the cyclic phase in F = 2 condensate is one more 
example of the systems, where spin-mixing is important. 
It is remarkable that, in this case, spin-mixing is happen- 
ing in the ground state and even in absence of rotation, 
applied magnetic fields and other external perturbations, 
as in F = 1 case. For the spin-mixing dynamics in i^ = 1 
condensate, see Ref. [23 • Note that recently, a new ex- 
perimental method for nondestructive study of internal 
degrees of freedom of spinor condensates was proposed 
in Ref. ^. 



FERROMAGNETIC STATE IN F 
CONDENSATE 



1 SPINOR 



In this Section we discuss thermal fluctuations in 
spinor F — 1 condensate. We concentrate on the ferro- 
magnetic state. The energy of the system is independent 
on the direction of the spin but it depends on the gra- 
dients of spin. The similar situation exists for the phase 
of the order parameter in scalar condensate (the energy 
is independent on its value, but depends on the gradi- 
ent), and therefore we can expect the similar behavior, 
namely, absence of long-range order in low dimensions. 
This result was obtained before for the large class of dis- 
crete and continuous spin models [l3|. We assume that 
in the ground state all the atoms occupy the only one hy- 
perfine state, ^_i = ^0 = 0, ^1 = y/n, and study small- 
amplitude oscillations induced by the temperature. In 
the ground state, spin is oriented along the z-direction. 



Jz — t, bx — 

of Sx and Sy : 



0. Perturbations of ^q lead to that 



5Sx 



SSy 



l2n 



(<5*o 



8%) 



(<5*o-^*S) 



(25) 



(26) 



The spectrum for the F =\ condensate was obtained 
in Refs. 113, 1"|. The frequency of the mode, responsible 
for fluctuations of ^0, is given by 



Till} = 



2m 



(27) 



Using Eqs. (25)- (27) we have calculated mean square 
fluctuations for Sx and 5*^ at large distances and found 
an absence of long-range order in ID and 2D cases. As an 
example, we show here the correlator in the 2D situation: 



[Sx{r) - Sx{0)y 



Td AT 



(28) 



where 



At — 



\/2mkBT 



(29) 



The derivation is similar to the case of cyclic phase in 
F — 2 condensate, presented in Section IV. The re- 
sult for Sy is the same. For ID case, the correlators 
behave as ~ r. This implies that in low-dimensional 
systems there is no long-range order in the direction of 
spin. However, the 2D system can be divided into blocks 
of a characteristic size L, £_s ^ L <^ Xt exp{Tci/T), 
with nearly the same direction of spin in each block 

{S,s ^ \ i^'kfi is the spin healing length). Different 
blocks are uncorrelated with each other. At temperatures 
of the order of 2Td/ \n{kBTd/ /3n), ^^ becomes compara- 
ble to At exp{Td/T). In this case, the size of each block 
becomes of the order of ^s . 

In the quasi-lZ? and 2D trapped condensates, which 
have finite sizes exceeding ^s, at low temperatures, one 
can expect the same orientation of spins throughout the 
system, and at higher temperatures, the cloud should 
consist of uncorrelated blocks. This is similar to the 
problem of fluctuating phase in quasi-lD and 2D scalar 
condensates 0, S 13 ■ Note that the case of ferromagnetic 
F — 2 condensate is similar to that of spin-1 system. 



VI. CONCLUSIONS 

We have studied the effect of thermal fluctuations on 
homogeneous infinite Bose-Einstein condensate with spin 
F — 2. We were interested in the cyclic state of this sys- 
tem, in which all the particles occupy three hyperfine 
states with mp — 0, ±2, and the energy depends on the 
relative phase x == 2So—S2~S-2 through the spin-mixing 
term. Using Bogoliubov-de Gennes equations, we have 
calculated mean square fluctuations of ^o, ^±2 and found 
absence of long-range order in one- and two-dimensional 
situations, but presence of this order for %. Then, we 
went beyond the small-amplitude oscillations approach 
and mapped our problem on the sine-Gordon model. A 



structure and energy of a single kink separating two topo- 
logical sectors with different values of x were found. In 
ID and 2D situations, thermal proliferation of kinks can 
lead to the roughening transition destroying a long-range 
order for x- One-dimensional infinite system is always in 
a rough phase, and finite ID system, as well as infinite 2D 
system, experience roughening transition at finite tem- 
peratures, whose values we have estimated. At the end, 
we have also discussed thermal fluctuations in ferromag- 
netic F = \ system and found an absence of long-range 
order in the direction of spin in low-dimensional situa- 
tions, as it can be expected from the Mermin-Wagner- 



Berezinskii theorem. We defined the typical size of each 
block, within which still there is an order in spin degrees 
of freedom. We expect that finite systems should be in 
the ordered phase at low temperatures. 
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VII. FIGURE CAPTIONS 

Fig. 1. The structure of a single kink in the cyclic 
state. Solid curve denotes x^ dashed one corresponds to 
So, and dotted curve shows 5*2 = S-2- 
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